Slow fracton modes in antiferromagnetic correlations of high-T c superconductors 
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In cuprate superconductors superconductivity develops as a unique crossover between the two 
, extremes characterizing these compounds, an antiferromagnetic (AF) Mott insulator phase, on 

one side, and a non-superconducting itinerant metal on the other. The hole doping tunes a 
cuprate inside such a broad spectrum of phases. Numerous inelastic neutron scattering stud- 
ies on underdoped samples shows most directly that the static AF order of the parent phase 
does not disappear by doping but transforms in such a way that at least fluctuating and local 
AF order persists, coexisting with superconductivity. Structurally, the model of-iptripe-type 
correlationsu accounts for the coexistence while their genuine dynamics might hideOu the most 
interesting physics of cuprates, including the pairing mechanism. Indeed, there is an impressive 
accumulation of knowledge in favour of AF correlations underlying the pairing. Some crucial 
properties of these coimelations have recently been pointed out: they should involve small ve- 
locity collective modestl of the wave vector identical to that of the static AF ordera iwhile their 
energy scaleO'u (< 50meV) should be rather low compared with any other involvedu. Here we 
show on phenomenological ground that the latter dynamics originates, irrespective of structural 
details of stripe order, from general geometrical constraints known to characterize the problem 
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^ 1 of disordered interpenetrated phases. In particular, we show that dynamics of the Euclidean- 

to-fractal crossover introduces a new (10-50 meV) energy scale, associated with slow (or almost 
localized) and well-defined fracton modes, and provides a consistent and natural interpretation 
for redistribution of spectral weight and other inelastic neutron scattering observations. We 
conclude therefore that the fracton excitations substantially contribute to pairing interactions 
'-^J in high-T c cuprates. 

o 

o : 

High-resolution studies of dissipation in weak link network of polycrystalline cuprates have recently identified pa 
fractal dissipative regime inside which a self-organized current-carrying medium is an object of fractal geometrya. 
Taking into account that the Josephson-coupled disordered medium represents a good description of the problem 
of phase coherence not only between the CuC>2 planes but also inside them (due to stripes and/or other forms of 
charge separation and clustering, filamentary fragmentation, etc.) the presence of fractal geometry in self-organized 
structures of CuC>2 planes seems rather plausible as well. In this work we investigate elastical modes in such a 
hypothetical (but still physically justified) two dimensional (2D) fractal lattice. In particular, we study AF spin 
0^ i fluctuations inside a model which allows a fractal geometry of the relevant AF cluster; the latter, being sample-sized 
and homogeneous in the parent compound, gets progressively more ramified by hole doping. (The liquid-crystal-like 
stripe phasesB represent of course a closely related visualization) . 

Generally, the dynamical features of such a problem are well-known: instead of extended modes (phonons or 
magnonsjj-of elastic quasi-continuum the vibrational modes in fractal lattices are localized excitations, known as 
fractonsBtJ. The fractality of real lattices depends on wave vector of the considered excitation so the fractoa, modes 
start to dominate only above some critical wave vector, i.e., above the frequency of phonon-to-fracton crossoveJ13. Also, 
there is a considerable accumulated knowledge about fractons specifically in classical diluted antiferromagnets both in 
their theoretical!! 2 ] and experimentaHliJ aspects. In these classical systems the magnetic component responsible for 
static AF order (e.g., Mn 2+ ion in MnF2 or in RbMnFa) is substituted at random by non-magnetic substituent (e.g., 
Zn 2+ or Mg 2+ , respectively). Inelastic neutron scattering studies on samples comprising substituent in concentrations 
approaching threshold concentration (above, which there would be no magnetic order) revealed dynamical features 
consistent with predictions for AF fractonst^lo. However, besides some doubts regarding proper assignment of specific 
experimental featurestJ there are still some uncertainties about the proper analytical form of the dynamical structure 
factor to be applied in modeling fractal antiferromagnets. 

Applying this knowledge to the problem of underdoped cuprates we first simplify our elaboration by limiting 
ourselves to the case of localized magnetic moments (spins associated to Cu 2+ ions) only. There is a consensus now 
that this is entirely appropriate for the spin wave phase of parent AF compounds but also equally inappropriate 
for the overdoped non-superconducting systems being compatible with itinerant electron description. The present 
work explores primarily the underdoped range. Localized description certainly makes sense for low hole concentration 
but becomes questionable close to optimal doping. There is however an increasing number of arguments^ pointing 
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out the localized nature of incommensurate magnetism even in optimally superconducting samples (oxygen doped 
La2Cu04 +J( ) so one could expect that the model we elaborate here may be applied at least to a reasonably wide range 
of underdoped concentrations. Now we define our 'building block' which, playing the role of a non-magnetic ion in 
diluted antifcrromagnets, are to be placed at random in the initially homogeneous and antiferromagnetic Cu02 plane. 
The fluctuating stripe correlations are certainly the most convincing interpretation for the incommensurate magnetism 
in superconducting cuprates and the local stripe order seems to be an inevitable element of cuprate physics. Our 
building blocks are therefore the regions with stripe order, disordered in sizes and, possibly, in other features. By 
'stripe order' we mean, in accordance with the localized magnetism of the present model, the originalj pattern of 
'the three spins wide AF stripes separated by non-magnetic antiphase domain boundaryld. The internal structure of 
the building block complicates somewhat the geometrical percolation scenario of the diluted antifcrromagnets; thus, 
we do not attempt to define the threshold hole concentration here. However, one expects that the main qualitative 
elements of an crossover to fracton dynamics, in increasing hole concentration, remain valid. 

Qualitatively, the scenario for underdoped cuprates would be the following: With no holes added the magnetic 
collective modes are high- velocity (Tiv — 650meVA, typically) regularly damped spin waves (Fig. la) characterized by 
dominantly 2D dynamics. The latter follows from the very high inplane superexchange (lOOmeV, typically) compared 
with the interplane one. Added holes introduce random obstacles for spin waves localizing them into fracton modes 
at short wave lengths (q > q c ), i.e., above the hole-concentration-dependent cross over energy Tiui c = Tivq c . This takes 
place when ramified AF cluster becomes self-similar (fractal) on all spatial scales smaller than l/<7 c j_subsequently, the 
crossover frequency progressively decreases as hole concentration (thus cluster fractality) increasesc3. 

Quantitatively, all of these dependencies must obey specific power-laws involving several statical and dynamical 
exponentsErEj (e.g., cocrelation length and dynamical one, fractal and AF fracton dimensions) which are, at least 
approximately, knownEj. The dynamical structure factor S(q, uj) has to include these power-laws and to assure known 
asymptotic behaviour characterizing the homogeneous (ui <C ui c ) and the fractal {v-p> u c ) regime. The problem of 
dynamical structure factor for AF fcaetons has been treated by many approaches!^ most of which are founded on 
the effective-medium considerations!!^!!!!-! Opt of several probably equivalent results for S(q,ui) in this work we use 
the form of Polatsek and Entin-WohlmanE£l as it explicitly demonstrates the localization features of fracton modes 
(Fig. lb). Although the latter result suffers £cpm general shortcomings of all effective medium approaches, why it 
cannot quantitatively be correct in all detailailJ, its general Green's function [background assures its basic usefulness. 
This form has been successfully used in studies of fractons in silica aerogelslij, and, more recently, in fused silicaEi 
Therefore, the form we propose to determine the scattering cross section of AF spin waves in cuprate CuC>2 planes of 
increasing fractality reads as 



S(q,u) 



v(u>)q 



r( w ) 



(ui - v{uo)q) 2 + T 2 (lu) (lj + v{uj)q) 2 + T 2 (u) 



(1) 



where v(uj) and r(w) are velocity of spin waves and damping (lifetime) function, respectively. The frequency depen- 
dence of these two quantities makes the present S(q,uj) different, together with its prefactor, from standard Lorentzian 
'shape functions '.-generally used in scattering problems and, particularly, in the inelastic neutron studies on diluted 
antiferromagnetscJ. The asymptotic forms of v(ui) and T(u>) are precisely knowntEt At low u> (ui < ui c ), v = vo, the 
magnon velocity, while at high lo (lo > lo c ), in fractal regime, it must obey v(ui) cx ui z (z is dynamical exponent). In 
the same limits the damping T(ui) crosses over from the phonon (magnon)-like Rayleigh law ( T(cu) oc u d+1 /uif, d is 
Euclidean dimension) to the localized Ioffe-Regel (r(w) oc ui) regimeca where both the dispersion and the wave vector 
description bepjmes ill-defined. For the specific frequency dependence of v and T we use the heuristic expressions of 
Courtens et.allij shown to work well in different fractal systemsEl Taking into account a predominantly 2D nature of 
magnetism in Cu02 planes the expressions are: 



v{u) = v [i + {u>/u c ) m ] z/m , r( u ) = r oUc ( w K) 3 [i + ( w Kn- 2/ra , (2) 

where m characterizes the sharpness of th^cKissover. Most of the features of the present model have a weak dependence 
on the choice of m and the ' traditional value m = 4 has beenjised. Also, in calculation of z = df/d a f, where 
dj is the fractal and d a f the AF fracton dimension, a well-acceptedt2l values df w 1.9 and d a f ~ 0.9 have been used. 
We note however that our results depend very little, at least qualitatively, on the fine tuning of exponents. Fig.l 
visualizes S(q,uj) of Eq.l taking into account real physical parameters of YBCO superconductor. 

Now we compare S(q,u) with the experimental results for imaginary part of generalized susceptibility y(q, w) on 
cuprates which determines the scattering cross section. These results are summarized in, e.g., Refs. |23|j2l[ Results 
at fixed q (w-scan) or at fixed u> (q-scan) depend very much on convolution with instrumental resolutionEj involving 
the integration of intrinsic S(q, ui) inside usually a sizable resolution ellipse. Hence, a better defined experimental 
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quantities are momentum- or energy-integrated .y ,(q, u>). The local (Brillouin zone averaged) susceptibility X2d( uj ) oc 
J x (<l,u)d 2 q is particularly well investigatecOcJ. In hole doped samples X2d( w ) snows a characteristic maximum 
(Fig. 2b, c) in the range of 20-50 meV while the parent compounds reveal approximately energy independent x'j^tiaEj, 
in agreement with S (q, uj) for spin waves. It is important to note that a sizable area of X2d( w ) P ea ks is positioned, 
in graphs involving absolute calibration, above the spine wave level. The latter reflects a substantial redistribution 
of spectral weight in frequency domain as the holes are added to the planes. The present model for S(q,u>) readily 
accounts for the redistribution. The 2D integration of our model's S(q,uj), Eq.l, gives: 



S 2 d(^)= I S(q,uj)d 2 q = 



where Q is Brillouin zone boundary and 



rAM^-lj+wln 



(i;Q-^) 2 + r 2 {vQ + uj) 2 + T 2 



r 2 uj 2 + r 2 



(3) 



, , vQ — lu vQ + UJ UJ 
A(ui) = arctg arctg h 2arctg-- (4) 

We note again that v and T are w-dependent (Eq.2). The result for S 2 d(w), which makes a central result of this 
paper, is shown in Fig. 2a for one particular choice of parameters together with the level which would characterize 
magnons propagating with velocity vq. Except behaviour just above the origin there is a full reproduction of ex- 
perimental curves. The physical reason for the redistribution in our model is clear immediately from Fig.l: For 
very small energies the magnetic modes start to propagate as magnons but soon, by approaching u> c , slow down and 
accumulate in the region around uj c taking away the spectral weight that magnons would have in the hole-free lattice 
up to high energies. We discover therefore that Tiuj c represents a natural new scale of energy characterizing the AF 
correlations in cuprates. S 2 d{w) of Eq.3 perfectly models most of the experimental results by choosing Tiui c from the 
interval 10 — 50meV. The present model attributes the experimental peaks in \2D ( an d also in x (<Z, w), see below) 
to dynamics of crossover quantified, in this work, by Eqs.2,3. One should note, however, that the main features are 
not determined by some special choice of spectral function S(q, uj), such as the one of Eq.l; instead, the features stem 
from universal principles of matching the two different regimes characterized by, different geometrical constraints. In 
our model tijjeije is no principal difference between the narrow 'resonant' peakcj, setting-in below T c , and the broad 
off-rcsonantElc3 peaks as they both reflect the dynamics of crossover. The experimental differences may come primarily 
from the involved parameters, such as the damping scale Tq. Indeed, it seems very plausible that the macroscopic 
superconductivity of the optimum-doped samples allows higher level of fractal self-organization accompanied with the 
increased value of T . One easily verifies that the increased r induces, in turn, a more pronounced peak in S^dO^). 

Now we analyze results for [fixed fixed q (w-scans) or at fixed uj (q-scans) in the light of our model, Fig. 3. In 
literature, there are many datac3 for x" with q fixed at AF wave vector, q = (0.5, 0.5) in reciprocal lattice units (rlu), 
measured on a variety of underdoped and optimally doped samples (no significant magnetic correlations 

has been observed in overdoped range). These results reveal a pronounced peak with the maximum positioned 
usually inside the interval 20-50 meV but with other details (e.g., temperature dependence, energy width) being 
dependent on the sample family and the extent of doping. In particular, the resonant peak remains as the only 
magnetic contribution to x m optimally doped samples (not yet observed in La2- a ;Sr 2 ;Cu04 family). In our model 
the region of Brillouin zone which belongs to the close vicinity of q = (0.5, 0.5) (say, for 0.5rlu < q < 0.6rlu, Fig. 3) is 
subject to an abrupt transformation of modes from a purely dispersive (spine wave) to almost localized (fracton) ones. 
A cj-scan taken nominally at q = (0.5,0.5) actually integrates the neighbouring region in which this transformation 
takes place. In the integration, the dispersive mode gives an (^-independent background contribution while the modes 
accumulating close to ui c contribute to the net w-dependence of x ((0.5, 0.5), uj). This is illustrated and further 
commented in Fig. 3. n[ — . 

Finally, we discuss the q-scans. The neutron data have shown that there is a characteristic linear relationshipclEZl 
between-.the q-space width of x (<Z> w ) and T c at constant uj . The slope of this dependence, having the dimension of 
velocityQ, has been usedu'0 as an evidence for the slow excitation, perhaps with a crucial contribution to formation of 
the superconducting state. In our model, a pronounced q-width of x (Qj w ) is a natural consequence of progressive 
slowing down and localization of magnetic modes close to ui c . The peak q-width depends on the value of u> c determined, 
in turn, by the hole concentration. Fig. 4 illustrates the increase of the q-width by decreasing uj c demonstrating that 
an approximately linear relationship is well-obeyed in a broad range of crossover frequencies. 

We conclude therefore that there is a high level of agreement between the present model and the experimental results 
for imaginary susceptibility and that the slow fracton modes could provide a basis for interpretation of the neutron 
studies. We also claim that these modes could provide a basis for understanding pairing in high-T c superconducto 
Generally, the fracton-mediated superconductivity has been treated theoretically in framework of quasi- crystals 
by solving an appropriate Eliashberg equation. On other hand, the coupling of carriers to 40 meV 'resonant spin 
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excitations' has recently been recognizedEl as strong enough to account for high values of T c in cuprates. The present 
work demonstrates that the fracton modes possibly underly these excitations and we believe that they could be a 
clue for 'uncovering the nature of slow moving charge objectstl, thus providing 'a crucial step for our understanding 
of high-T c superconductivity '0. 
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FIG. 1: Dynamical structure factor S(q,ui) for inelastic neutron scattering by cuprates' AF correlations of this model, Eqs.1,2. 
The two asymptotic dynamics is shown: dispersive spin waves (shown in a)) and localized fractons (shown in b)). Lattice 
parameters of YBCO superconductor have been assumed. Known velocity of spin waves vo has been used as a parameter. 
Some reduction of vo by hole doping has been both observed and expected. The form of Eq.2 has been used for damping; 
however, in a) the damping was somewhat additionally suppressed in order to visualize the spin waves. The main physical 
parameter is the crossover frequency lj c . It decreases by hole doping and defines the energy scale of localizing fracton modes. 
The feature of localization is reflected by vanishing wave vector dependence for lo > lo c , where the mode extension becomes 
parallel to q-axis. 
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FIG. 2: Integrated (local) dynamical structure factor S2d{u) of Eqs.3,4, normalized to 6*215 of AF magnons characterizing the 
Euclidean parent phase, is shown in a) (full black line). A peak develops as a result of accumulating modes with frequencies 
close to the crossover (lu c ) one. Numerical integration of Eq.l has also been performed (open circles). Disagreement for small 
energies represent an artifact of numerical integration of very sharp functions. Experimental results for local susceptibilities of 
doped superconducting cuprates are shown in b), Ref.23, and c), Ref.24,26 on the same energy scale. Thick gray lines in all 
panels represent the level associated to propagating spin waves. 
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FIG. 3: tu-dcpcndence of S(q, u>), Eqs.1,2, for five fixed wave vectors close to the AF zone center at qAF = 0.5. The numbers near 
the curves designate the parametric values of q in reciprocal lattice units 2n/a oc 1.63A - . A crossover to a weak g-dependence 
(i.e., localization) obviously takes place in this interval of q. Thick gray line is the average of curves. Depending on parameters, 
there are sharp propagating magnon peaks (not shown) immediately above qAF (0.5 < q < 0.55, typically). Inset: experimental 
result on underdoped superconductor, Ref.26 
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FIG. 4: Dependence of the relative g-width of S(q, cu) on doping-dependent parameter ui c . Inset shows a set of u;-scans of S(q, u>), 
Eqs.1,2, for q fixed at q — 0.6 differing only in the designated choice for u> c . The amplitude of 5(0.6, a;), taken at Tiui = 25meV , 
has been chosen as a relative measure for the g-width. One notes a broad range of lo c characterized by approximately linear 
behaviour of g-width. Full line is a guide for the eye. 



